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CN ■ ABSTRACT 

' This paper presents a first observational investigation of the faint Of?p star NGC 1624-2, 

yielding important new constraints on its spectral and physical characteristics, rotation, mag- 
£vq ■ netic field strength, X-ray emission and magnetospheric properties. Modeling the spectrum 

t— I ' and spectral energy distribution, we conclude that NGC 1624-2 is a main sequence star of 

. mass M 30 M , and infer an effective temperature of 35 +2 kK and logg = 4.0 + 0.2. Based 

. ■ on an extensive time series of optical spectral observations we report significant variability of 

\ a large number of spectral lines, and infer a unique period of 157.99 + 0.94 d which we inter- 

■ pret as the rotational period of the star. We report the detection of a very strong -5.35 + 0.5 kG 

- longitudinal magnetic field (B-), coupled with probable Zeeman splitting of Stokes / pro- 
files of metal lines confirming a surface field modulus (B) of 14 ± 1 kG, consistent with a 
surface dipole of polar strength > 20 kG. This is the largest magnetic field ever detected in an 
O-type star, and the first report of Zeeman splitting of Stokes / profiles in such an object. We 
also report the detection of reversed Stokes V profiles associated with weak, high-excitation 
emission lines of O m, which we propose may form in the close magnetosphere of the star. 
We analyze archival Chandra ACIS-I X-ray data, inferring a very hard spectrum with an X-ray 
efficiency log L x /Lboi = -6.4, a factor of 4 larger than the canonical value for O-type stars and 
comparable to that of the young magnetic O-type star 9 l Ori C and other Of?p stars. Finally, 
we examine the probable magnetospheric properties of the star, reporting in particular very 
strong magnetic confinement of the stellar wind, with rj* =* 1.5 x 10 4 , and a very large Alfven 
radius, /? A if = 11.4/?*. 
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* Based on spectropolarimetric observations obtained at the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the National Re- 
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1 INTRODUCTION 



The detecti on (e.g. Grunhut et al.l |2009|) , empirical characteri- 
zation (e.g. I Wade et . al.l 1201 lh and theoretical modelling (e.g. 
ISundqvist et aTT2012l) of a growing sample of magnetic O-type 
stars is leading to a new, refined picture of the scope and impact 
of magnetic fields in high-mass stars. 

O-type stars are unique laboratories for investigating the 
physics of stellar magnetism. Magnetic fields have clear influ- 
ence on their rotation rates (rotation periods of most detected 
magnetic O stars are significantly longer than those of non- 
magnetic O stars of similar spectral types; e.g. lUd-Doula et al.l 
l2009MHowarth et al.l200l;lMartins et alfefJlOl) . Evolutionary mod- 



Table 1. Log of spectroscopic observations used to determine the pe- 
riod. B refers to the spectral range that includes He i /1447 1 , He n /14542, 
He n ,14686, and H/3; V to the range that includes He i /15876; and R to the 
range that includes Ha. J stands for the low-resolution observations from 



els dMaeder & Mevnetl l200j |2004> and, recently, observations 
(Briquet et al., MNRAS, submitted) of massive stars suggest that 
the internal rotation profile is strongly modified by the presence of 
a magnetic field, enforcing essentially solid-body rotation through- 
out the bulk of the outer radiative zone. Magnetic fields have clear 
and fundamental effects on the structure, dyn amics and heating of 
the powerful radiative wind s of O stars (e.g. lud-Doula & Owoclq 
l2002MSundqvist et al.l2012h . The lives of magnetic O-type stars are 
therefore expected to differ significantly from those of their non- 
magnetic brethren. 

The subject of the present paper, NGC 1624-20, is an Of?p 
star dWalborn et alj20ld) and the main ionizing source of the open 
cluster NGC 1624 (the H n region S212). It is also one of only 
eight O-type stars in which magnetic fields have been dete cted with 
confid ence. The classification Of?p was introduced by Walborn 
(1972) to describe spectra of early O-type stars exhibiting the pres- 
ence of C in /14650 emission with a strength comparable to the 
neighbouring N m lines. Well-studied Of?p stars are now known 
to exhibit periodic spectral variations (in Balmer, He i, C ni and 
Si m lines), narrow P Cygni or emission components in the Balmer 
lines and He i lines, an d UV wind l ines w eaker than those of typ- 
ical Of supergiants (see lNaze et al.l J2010l) and references therein). 
With our report of a detection of a magnetic field in NGC 1624- 
2, magnetic fields have now been firmly detected in all of the 5 
known Galactic members o f this class - HP 191612: iDonati et al.l 
J2006h . IWadeetal.1 (1201 ll). HP 108: iMartins et al.l feoioh ; HP 
148937: IHubrig et al] d2008h . IWade et all d2012l) : CPP -28 2561: 
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- prompting the inference that there is a direct physical relationship 
between the magnetic field and the Of ?p c haracterist i cs. 

According to the recent analysis by Ijose et al"] d201ll) , NGC 
1624 is a young open cluster located significantly above the Galac- 
tic plane. Their analysis yields a heliocentric distance of 6.0 + 
0.8 kpc (in agreem ent with that found in the pioneering study of 
lMoffatetalJll979l) and an age of no greater then ~ 4 Myr. In addi- 
tion to NGC 1624-2, which is by far the brightest cluster member, 
3 other apparently bright optical sources are located within 2 ar- 
cmin of the cluster centre. These sources are reported to be an early 



the University of Hawaii, as well as on observations obtained using the Nar- 
val spectropolarimeter at the Observatoire du Pic du Midi (France), which is 
operated by the INSU. The spectroscopic data were gathered with five facil- 
ities: the 9.2 m Hobby-Eberly Telescope at McDonald Observatory (MDO), 
the 3.5 m Telescope at Calar Alto Observatory (CAHA), the 1.5 m Tele- 
scope at the Observatorio de Sierra Nevada (OSN), the 4.2 m William Her- 
schel Telescope at the Observatorio del Roque de los Muchachos (ORM), 
and the 2 m Himalayan Chandra Telescope at Indian Astronomical Obser- 
vatory (IAO). 

t E-mail: wade-g@rmc.ca 

1 According to the numbering system of lMoffat et alj ll979l) . 



Jose et l|), G for the intermediate-resolution GOSSS observations 

iMaiz Apellaniz et al]|201ll), and N fo r the high-resolution NoMaDS ob- 
servations iMafz^peliarriz et al.l201 ll) . 
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B-type main sequence star, and two probable F giants (which are 
probably not physically associated with the cluster). 

The current paper provides a first observational characteriza- 
tion of NGC 1624-2, the faintest of the known Galactic Of?p stars 
(V = 11.8). In particular we report the physical parameters of the 
star, its rotational period, the detection of its magnetic field, and 
a preliminary characterization of its magnetic and magneto spheric 
characteristics. In Sect. 2 we describe the spectroscopic and spec- 
tropolarimetric observations upon which we base our results. In 
Sect. 3 we provide a short overview of the spectral properties of 
the star. In Sect. 4 we derive the physical properties of the star and 
its wind. In Sect. 5 we determine the spectral variation period of 
the star. In Sect. 6 we describe the magnetic field diagnosis and our 
constraints on the surface magnetic field. In Sect. 7 we characterize 
its X-ray properties based on archival data. In Sect. 8 we describe 
the probable magnetospheric properties of the star. In Sect. 9 we 
discuss our results and summarize our conclusions. 



2 OBSERVATIONS 

2.1 Spectroscopic observations 

The optical spectroscopic observations for the period determina- 
tion were obtained within the context of three different projects. 
High resolution (R = 30000) spectroscopy was obtained with 
the NoMaPS project dMafz Apellaniz et all 1201 ll) using the 9.2 
m Hobby-Eberly Telescope (HET). Intermediate resolution (R = 
2 500-4 000) spectroscopy w as obtained within the GOSSS project 
dMafz Apellaniz et alj|201 lb using the 1.5 m telescope of the Ob- 
servatorio de Sierra Nevada (OSN), the 3.5 m telescope of the Calar 
Alto Observatory (CAHA), and the 4.2 m William Herschel Tele- 
scope (WHT). Finally, low resol ution (R = 1 500) spectroscopy 
was obtained by I Jose et alj d201lh using the 2 m Himalayan Chan- 
dra Telescope. 

The 33 NoMaPS spectra were gathered using the high reso- 
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Figure 1. Two NoMaDS spectra of NGC 1624-2 (top, obtained on 22 Aug 2011, high state, phase 0.92; middle, obtained on 11 Nov 2011, low state, phase 
0.43). The bottom spectrum represents the difference (high minus low). The spectra have been convolved to a resolving power R = 10 000 for display purposes. 
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Figure 2. Of?p-diagnostic emission lines of C in and N m near maximum 
emission at the full ESPaDOnS resolving power (R = 65 000). Note the re- 
markable composite emission profiles with broad and narrow components, 
not seen in previously investigated members of the Of ?p category. 



lution fiber-fed echelle spectrograph (HRS) at the HET. The spec- 
tra were obtained with two different instrument setups (600g4739K 
and 600g6302) allowing us to cover the entire wavelength range of 
~3800-7300 A. The fiber sdiamter was 2". We obtained two ex- 
posures of 3600 seconds in the blue-setup (~3800-5700 A) and 
two of 1800 seconds for the red-setup (-5400-7300 A). The data 
were reduced following the IRAF/package-based procedure sug- 
gested on the HET we bsitjl. Cosmic rays were removed using the 
L.A. Cosmic package Jvan Dokkumll200lh . The spectra were flux 
normalized as well as wavelength calibrated. The typical resulting 
signal-to-noise ratio (SNR, per resolution element) is between 150 
and 250, depending on the wavelength range. For NGC 1624-2, ad- 
ditional spectra were obtained in order to remove the nebular con- 
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tribution. For these observations the fiber was positioned about 14" 
NE and SW of NGC 1624-2 to avoid stellar contamination, and the 
spectra were reduced as described above and subtracted from the 
stellar observations. 

The 19 GOSSS spectra were obtained with various spectro- 
graphs fed by the OSN, CAHA and WHT telescopes as part of 
the larger Galactic O-Star Spectroscopic Survey. The data were re- 
du ced using a pipeli ne specifically written for the project, described 
bv lSotaet all 1201 lh . 

A log of the 64 spectroscopic observations is presented in Ta- 
ble!]] where we specify which regions were observed on each date. 



2.2 Spectropolarimetric observations 

Five high resolving power (R 65 000) spectropolarimetric 
(Stokes / and V) observations of NGC 1624-2 were collected with 
ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) be- 
tween Feb 1 and 9 2012. An additional observation was obtained 
with Narval (essentially a twin of ESPaDOnS) on the Bernard Lyot 
telescope (TBL) on Mar 24 2012. All of the spectropolarimetric 
observations were obtained within the context of the Magnetism in 
Massive Stars (MiMeS) Large Programs. Each spectropolarimetric 
sequence consisted of four individual subexposures, each of 600 s 
duration (for ESPaDOnS) or 1 200 s duration (for Narval), taken 
in different configurations of the polarimeter retarders. From each 
set of four subexposures we derived Stokes I and St okes V spec- 
tra fol lowing the double-ratio procedure described by bonati et al.l 
jl997t) . ensuring in particular that all spurious signatures were re- 
moved to first order. Null polarization spectra (labeled AO were cal- 
culated by combining the four subexposures in such a way that po- 
larization cancels out, allowing us to verify that no spurious sig- 
nals are present in the data (see iDonati et al.lll997l for more de- 
tails on the definition of N). All frames were processed using th e 
automated reduction package Libre ESpRIT l lDonati et al.lll997T ). 
The peak SNRs per 2.6 kms~' velocity bin in the reduced spec- 
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Figure 3. Left panel: Fit of the observed He lines (black solid line) by three CMFGEN models with r c ff = 33 000 (green long-dashed line), 35 000 (red 
short-dashed line) and 37 000 K (blue dotted line). Right panel: Best SED (black solid curve) from the CHORIZOS run with fixed T eff = 35 000 K. The 
triangles with error bars are the photometric magnitudes used in the fit (UBVIJHK), the circle represents the measurement that was unused (R), and the stars 
are the synthetic magnitudes corresponding to the best SED. 



tra range from 86-140, with the variation due principally to seeing 
conditions. 



3 THE OPTICAL SPECTRUM OF NGC 1624-2 

Fig. Q] shows the blue-green spectra of NGC 1624-2 from the 
HET/NoMaDS data (convolved with a Gaussian to a uniform R = 
10000 for classification purposes) near maximum and minimum 
emission phases, together with their difference. The apparent spec- 
tral types from the He ionization ratios are 06.5 and 08, respec- 
tively, although both are affected by particular spectral features: at 
maximum, longward emission filling in He i ,44471 is clearly seen, 
while at minimum, the He i and He n profiles are very different due 
to the stronger Stark wings in the latter, and indeed, the equivalent 
width ratio ~ 1 corresponds to a type near 07. 

Fig.|2]displays the N in /14640, C m /14650 complex near max- 
imum emission at the full ESPaDOnS resolving power of 65 000. 
NGC 1624-2 exhibits the most extreme Of?p spectrum known to 
date, with the C ra emission lines stronger than those of N m, 
and stellar Balmer emission much stronger than in any other. The 
N m and C ni features (as well as, we suspect, He n /14686) dis- 
play remarkable composite emission profiles with broad and nar- 
row components, not seen in previously investigated members of 
the Of?p category. We hypothesize that the broad components ar e 
photospheric as in normal Of or Ofc stars dWalborn et al.ll2010l) . 
while the narrow components typical of Of?p spectra are magne- 
tospheric. In agreement with this interpretation, it is seen in Fig.Q] 
that only the narrow components vary with phase. However, the 
narrow C m emission does not disappear at minimum, so unlike 
most other members of the category NGC 1624-2 remains Of?p at 
both extreme phases. Moreover, the apparently greater intensity of 



C m as compared to N m is largely due to the broad component, so 
that the intensities of both species in the narrow system are more 
comparable, similarly to other Of?p spectra. 

This is the first case combining Ofc and Of?p emission 
profiles, for which we introduce the new notation Of?cp; thus 
the full classifications of these spectra are now 06.5f?cp and 
Q8f?cp at max i mum and minimum, respectively. As discussed by 
IWalborn et al.l | |2010|) . the Ofc phenomenon in normal Galactic 
spectra is strongly peaked toward type 05. However, the optical 
component of the high-mas s X-ray binary HP 153 919 (2U 1700- 
37) is 06.5 Iafc (see, e.g., IConti & Cowievlll975l, A. Sota et al. 
in prep.), as i s the primary of the LMC double-lined spectroscopic 
binary R139 dTavlor et alj|201 if) . Thus, the Ofc phenomenon also 
occurs at later spectral types in some anomalous objects. 



4 STELLAR PHYSICAL AND WIND PROPERTIES 

The atmospheric parameters of NGC 1624-2 have been derived by 
means of spectroscopic analysis with the code CMFGEN, fitting 
the minimum-emission state spectrum (the average of NoMaDS 
spectra acquired on 17/18 Nov 2011) . A fu ll description of CM- 
FGEN is provided bv lHillier & Milled d 19981) . In brief, it computes 
non-local thermodynamic equilibrium (non-LTE) atmosphere mod- 
els including winds and line-blanketing. The statistical equilibrium 
and radiative transfer equations are solved in the comoving frame 
through an iterative scheme to compute the level populations and 
radiation field. The temperature structure is set by the condition of 
radiative equilibrium. Our models include the following elements; 
H He, C, N, O, Ne, Si , S, Ar, Ca, Fe, Ni. The solar abundances 
of lGrevesse et al.l d2007l) are used unless stated differently. The ve- 
locity structure is specified at the beginning of the computation. It 



NGC 1624-2: A slowly rotating, X-ray luminous Of?cp star with an extraordinarily strong magnetic field 



is usually a combination of a pre-computed photospheric structure 
connected smoothly to a so-called /? velocity law (v oc (1 - r/R,f 
where r is the radial coordinate and R, the stellar radius). This 
structure is subsequently iterated to satisfy the momentum equilib- 
rium equation in the inner atmosphere. Once the atmosphere model 
is converged, a formal solution of the radiative transfer equation is 
performed to yield the emergent spectrum. A microturbulent ve- 
locity varying with height from 10 to 200 km s -1 (beginning in the 
photosphere and continuing to 100 R B ) is used. 

Prior to a comparison between synthetic and observed spectra, 
it is necessary to determine the projected rotational velocity and 
macroturbulent bro adening. We u s ed the Fourier transform method 
for the former (e.g. |jankovlll995l ; ISimon-Dfaz et al.ll2006l) . focus- 
ing on the C iv /15801 line. The Fourier transform did not show 
any zero before reaching the noise level, implying an upper limit 
on vsin! of ~ 15 km s" 1 . Adopting this value, we convolved our 
synthetic spectra with a rotational profile and a Gaussian profile 
aimed at reproducing isotropic macroturbulence. We found that a 
macroturbulent velocity of about 25 km s" 1 gave the best fit to 
the observed He lines. We note based on the results of Sect. [6] 
that this additional broadening is likely to be magnetic in origin, 
and so isotropic macroturbulence will provide only an approximate 
reproduction of its effects. The combined line broadening (rota- 
tional, turbulent, magnetic) of NGC 1624-2 is compa rable to that of 
the ve ry sharp-lined magnetic 09IV star HD 57682 dGrunhut et al.l 
l2009l and Grunhut et al., MNRAS, submitted) . 

The He i/He 11 ionization balance was used to derive T eS . As 
for the other Of?p stars, we relied on the spectrum closest to the 
minimum-emission state. It has the least contaminated He 1 lines. 
Based on the absence of significant variability of the majority of 
the He 11 lines in our spectra, we assume these lines are principally 
photospheric. Fig.[3](left panel) shows the He lines (7 lines of He 1, 
3 lines of He n; He 11 /14686 was not used, as it is strongly variable) 
used in our analysis together with three selected models at 33 kK, 
35 kK and 37 kK. Overall, the 35 kK model provides the best repre- 
sentation of the observed profile. The He 11 lines are either too weak 
or too strong if T c ff = 37 kK or 33 kK respectively. The He 1 lines 
also appear too strong at low T c f[. However, we note that the exact 
shape of the He 1 lines is uncertain (due to residual contamination 
by unmodelled contributions from the non-spherical component of 
the wind) and consequently, a safe uncertainty on T c g is 2000 to 
3000 K. 

The surface gravity of NGC 1624-2 was derived from the 
shape of the Balmer-line wings, as is usually done for O stars. We 
found that a value of log g = 4.0 was the best compromise to repro- 
duce all lines. A slightly higher value (log g = 4.2) is indicated by 
H/3 and He compared to Hy. We can safely exclude any value be- 
low 3.7, the broadening being insufficient to reproduce the observed 
profiles in that case. We thus adopt a value of log g = 4.0 ± 0.2. 

We also investigated t he nitrogen content of NGC 1624-2 by 
means of its N m lines (e.g. lMartins et alj|2012l) . We used the N m 
A4195, N m A4510-4525 and N m ,14535 features, as all other ni- 
trogen lines are either in emission, absent, or of very low SNR (be- 
low 10). Assuming the adopted microturbulence described above, 
a v alue of N/H=(1.0 ± .4) x 10~ 4 (by number) gives a satisfactory 
fit. lMartins et all J2012h tested the effect of a change of the micro- 
turbulent velocity in the atmosphere model computation (as well as 
the output synthetic spectrum) and found that a 5 kms~' increase 
in the microturbulence resulted in a 10-20% change in the derived 
abundance. This effect, while not negligible, is smaller than the ef- 
fect of r cff uncertainties. Ultimately, while the exact value of the N 
abundance is uncertain due to the turbulent/magnetic and tempera- 



Table 2. Results of the CHORIZOS modelling of the spectral energy distri- 
bution. The first column gives the results for the run in which the luminosity 
class (LC) was allowed to vary, but in which the effective temperature (7" e (f ) 
was fixed to 4.8, i.e. close to luminosity class 5. The second column gives 
the results for the ran in which the effective temperature was allowed to 
vary, but in which the luminosity class (LC) was fixed. 



Quantity 


Free LC 


Free T e g 


y 2 

'Wed 


2.03 


0.86 


7W (K) 


35 000 


31 000 ±2800 


logg(cgs) 


3.891 ±0.140 


3.957 ± 0.022 


log L (solar) 


5.125 ±0.171 


4.632 ± 0.214 


£(4405 - 5495) (mag) 


0.802 ± 0.015 


0.794 ± 0.018 


E(B - V) (mag) 


0.802 ± 0.014 


0.792 ± 0.016 


K5495 


3.752 ± 0.092 


3.736 ± 0.096 


A v (mag) 


3.040 ±0.031 


2.991 ±0.046 


logd(pc) 


3.712 ± 0.094 


3.537 ± 0.059 



Table 3. Summary of adopted stellar, wind, magnetic and magnetospheric 
properties of NGC 1624-2. The effective temperature T c ff, surface gravity 
logg, radius i?*, luminosity L+ and mass M+, as well as the reddening 
E(B - V), extinction Ry and distance log d were determined based on mod- 
elling of observations described in Sect. [4] The upper limit on the projected 
rotational velocity - determined from the inferred radius and period assum- 
ing rigid rotation - is consistent with the negligible additional line broad- 
ening required to fit the magnetically-split profiles of the C iv line profiles 
(Sect.|6). The wind mass-loss rate M and terminal velocity v m are calculated 
based on theoretical considerations in Sect. [8] The magnetic field dipole 
surface strength was estimated based on measurement and modelling of 
Zeeman effect, described in Sect. [6] Finally, the wind magnetic confinement 
parameter rj+, the Alfven radius i?Alf> the rotation parameter W and the Ke- 
pler radius /?Kep are calculated based on measured and infen'ed quantities in 
Sect. [8] 



Spectral type 


06.5f?cp-08f?cp 


7W (K) 


35 000 ± 2000 


log g (cgs) 


4.0 ± 0.2 


R*(Ro) 


10 ±3 


vsin 1 (kms -1 ) 


< 3 


log(L*/L ) 


5.10 ±0.2 


MT C (M e ) 


34 ±31 


M e ™ 1 (M ) 


28^ 


logMOVIoyr" 1 ) 


-6.8 


v OT (kms -1 ) 


2875 


E(B - V) (mag) 


0.802 ± 0.02 




3.74 ±0.1 


\ogd (pc) 


3.712 ±0.1 


Bd (kG) 


~ 20 




1.5 X 10 4 


R Alf (R*) 


11.4 


w 


4 X 10~ 3 


R K ep(«.) 


> 40 


T W m (Myr) 


0.24 



ture uncertainties, we can conclude that the star is not strongly N 
enriched ([N/H]<0.3). 

In the absence of UV spectroscopy, and given the peculiar 
shape of the emission in wind sensitive optical lines (He 11 /14686, 
Ha), it is not possible to derive the mass loss rate of NGC 1624-2 
based on the observations. First, the terminal velocity is not known 
(it is usually derived from the blueward extension of UV P-Cygni 
profiles). Secondly, the optical emission lines are very narrow, as 
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in the other Of?p lines, and are probably the result of a complex 
(magnetically-con fined) wind geometry t hat cannot be reproduced 
by our ID models llSundqvist et al.ll2012h . 

In order to derive the distance and extinction of NGC 1624- 
2 we used the Bayes ian (spectro)photometric code CHORIZOS 
dMai'z Apelliiuzll2004h . In the latest CHORIZOS version, the user 
can select distance to be an independent pa rameter by applying at- 
mosp here models (TLUSTY for OB stars; lLanz & Hubenvll200l 
|2007|) calibrated in luminosity with the help of Geneva stellar evo- 
lutionary tracks (excluding rotational effects). In such models the 
intrinsic parameters are effective temperature (T^ff) an d photomet- 
ric luminosity class (LC, defined in a similar way to its spectro- 
scopic counterpart, with 0.0 corresponding to hypergiants and 5.5 
to ZAMS). The extrinsic parameters are reddening (£(4405-5495), 
the monochromatic equivalent of E(B - V)), extinction law (R5495, 
the monochromatic equivalent of R v ), and logarithmic distance 
(log d). CHORIZOS uses the family of extinction laws of Mafz 
Apellaniz et al. (20 12, in preparation), which improve upon the 
ICardelli etalH l989) laws by using spline interpolation in the opti- 
cal range instead of a seventh-degree polynomial and by correcting 
the extinction in the U band. 

As inpu t photo metry, we used the UBVI magnitudes 
of Ijose et all d201lh and the JHK magnitudes o f 2MASS 



Skrutskie et al.l2006 ). We discarded the R magnitude of ljose et al .1 



201 11) due to the influence of the Ha emission on the stellar SED, 



though we later checked that its inclusion would not have greatly 
changed our results. We did two different runs: in the first one, we 
fixed T e ff to be 35 000 K (the value derived from spectroscopy) and 
allowed LC to vary in such a way that 3.7 < logg < 4.3 (slightly 
larger than the range derived from spectroscopy). In the second run, 
we left r c ff as a free parameter and fixed LC to be 4.808 (i.e. close to 
the spectroscopic luminosity class V), since that value corresponds 
exactly to logg = 4.00 (cgs) for r cff = 35 000 K and logg is ap- 
proximately constant for a fixed LC for O stars near LC = 5.0. In 
the two runs we left the three extrinsic parameters (reddening, ex- 
tinction law, and logarithmic distances) free. 

The results of the two CHORIZOS runs are shown in Table [2] 
and the mean SED from the first run is shown in Fig. [3] The two 
runs show good values of the reduced x 1 -, indicating that the pho- 
tometry is consistent with the input SED models and extinction 
laws. All of the results are consistent (within 2cr) between the two 
runs, with the values for £(4405 - 5495), E(B - V), R 5495 , and 
A\ practically indistinguishable. The most significant differences 
between the two runs are in T c ff , log L, and log d: the first run indi- 
cates a hotter, less luminous, and closer object in comparison with 
the second run. The second run has a better xL. but that can be 
ascribed to a random fluctuation in the input photometry (the ob- 
tained T eB is within 2cr of the spectroscopic value) or a small effect 
of the circumstellar material in the SEDj. The distance obtained in 
the first run (5.2 ± 1.1 kpc) is in better agreement with the values 
from the literature (~ 6 kpc). The values for the extinction indi- 
cate that Ay is close to 3 magnitudes and that the value of ^5495 is 
slightly larger than the mean Galactic value of 3.1-3.2, something 
that is typical for objects in H11 regions. 

Using the CHORIZOS-derived distance, reddening and ex- 



3 We did equivalent CHORIZOS runs with the lCardelli et alj Jl989h ex- 
tinction laws and we found the same effect as in our derivation of the new 
Mafz Apellaniz extinction laws: for the fixed 7" c ff case^j is larger and for 
the free T^g case the measured value is further offset from the spectroscopic 
one. 
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Figure 4. Illustration of the variation of the He 1 /15876 line in the NoMaDS 
spectra. Numbers associated with each spectrum indicates date in the format 
yymmdd, and the phase computed according to Eq. (1). 



tinction, the best CMFGEN fit to the SED corresponds to 
logL/L Q = 5.08. With the quoted uncertainty on d, the uncer- 
tainty on logL/L is about 0.15 dex. Therefore the CMFGEN 
(5.08 ± 0.15) and CHORIZOS (5.13 ± 0.17) values of log L/L a are 
in excellent agreement. 

Based on the derived temperature (T c g = 35 000 ± 2000 K 
and luminosity (averaging the CMFGEN and CHORIZOS results, 
log L,/L = 5.1±0.2), we estimate the stellar radius R, = 10±3 R B . 
The mass of NGC 1624-2 derived from the spectroscopic analysis 
(i.e. from logg and the derived radius) is M, = 34 + 31 M Q . The 
uncertainty is very large because the error on log g is nearly a fac- 
tor of 1.6 (0.2 dex). This "spectroscopic" mass can be compared 
to the "evolutionary" mass derived simply from the (r eff ,L.) po- 
sition of NGC 1624-2 on the HR d iagram. Using the evolutionary 
tracks of iMevnet & Maeded ( feOOSl) . we find M evol = 28^ M Q . This 
estimate implicitly assumes that the evolutionary tracks used can 
represent the evolution of the star. Given the peculiar properties of 
NGC 1624-2 (in particular its very strong magnetic field, very slow 
rotation and probably modified mass loss) this is not completely 
obvious, and thus the evolutionary mass should be regarded as only 
indicative. 

The adopted parameters for NGC 1624-2 are summarized in 
Tabled 
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Figure 5. Left panel: Periodogram of the He n A4686 equivalent width measurements. Right panel: Equivalent widths of the He n /14686 line of NGC 1624-2, 
phased according to Eq. (1). The ESPaDOnS observations of this star were acquired at phase ~ 0.96, while the Narval observation was acquired at phase 
~ 0.28. 



Table 4. Equivalent width of He n /14686) (in A) in emission as a function 
of Julian Date. The final column indicates the origin of the spectrum, as in 
Tabled 



JD 


EW 


0"EW 


Type 


(-2450 000) 


(A) 


(A) 




3 985.8 


0.578 


0.035 


J 


3987.0 


0.502 


0.043 


J 


4127.6 


0.549 


0.026 


J 


4754.6 


0.673 


0.009 


G 


5 133.6 


0.665 


0.004 


G 


5 134.6 


0.723 


0.006 


G 


5 135.5 


0.658 


0.007 


G 


5 137.5 


0.766 


0.013 


G 


5 139.7 


0.732 


0.013 


G 


5 153.7 


0.898 


0.024 


J 


5 160.5 


0.950 


0.010 


G 


5 795.9 


1.036 


0.008 


N 


5 836.8 


0.932 


0.009 


N 


5 837.8 


0.986 


0.007 


N 


5 876.7 


0.531 


0.011 


N 


5 883.7 


0.514 


0.007 


N 


5 884.7 


0.530 


0.007 


N 


5 969.5 


1.041 


0.116 


G 


5 973.5 


1.034 


0.058 


G 


5 998.6 


0.918 


0.007 


N 



5 SPECTROSCOPIC VARIABILITY AND PERIOD 

The first indication of variability in the spectrum of NGC 1624- 
2 resulted from a comparison of the different GOSS S observations 
obtain ed in 2008 and 2009 (s ee Table [Tj>, reported bv l Walbqrn et al.l 
d2010h and lSota et all d201ll) . Indeed JWalborn et all d2010h hinted 
at the possibility of variability on a time scale of days and suggested 
the need for further observations. 

In order to determine the period of NGC 1624-2, we ana- 
lyzed different optical lines and line ratios. We settled on four 



lines: Ho-, Hei /15876, H/3, and He n /14686; and one line ratio, 
He ii /14542/Hei /14471. Of these six lines, three remain in emis- 
sion in all of our observations (Ho-, H/3, and Hen /14686), two 
remain in absorption (Hen /14542 and Hei /14471) and the last 
one (Hei /15876) changes from emission to a complex absorp- 
tion/emission profile. The variation of the He i /15876 line is il- 
lustrated in Fig. [4] 

We measured the emission equivalent width (EW; i.e. the in- 
tegral of <Mx(flux-continuum)/continuum)) of all the lines at every 
epoch and we analyzed the variations of the first four and the line 
rat io of the last two by mean s of Scargle periodograms as described 
bv lHorne & Baliuna s U1986h . We examined periods in the range 60- 
400 days, using 10 5 periods uniformly distributed in log (period). In 
all five cases the periodograms show maxima between 156 and 159 
days (Tabled- The periodogram for He n /14686 is shown in Fig. [5] 
(left panel). If we phase the EW measurements of this line with the 
adopted period (158 d, discussed below, and corresponding to the 
dominant peak of the peoriodogram), we obtained the variations 
shown in Fig.[5](right panel), where phase has been set for 9 Feb 
2012 at noon UT (i.e. HJD 2455967.0; these same measurements 
are summarized in Table |4). The behaviour is the same for all the 
five quantities measured, with maxima (in either the EW or the line 
ratio) around phase 0.10-0.25 and minima around phases 0.6-0.75. 
In addition, variations at the 0.01 magnitude level were detected in 
an /?-band monitoring program of NGC 1624-2 with the 0.6 m tele- 
scope at the Esteve Duran Observatory that are consistent with the 
changes in EW(Ha) detected with spectroscopy. 

To obtain a "best" period we tried two different approaches. 
First, we combined the five measured periods weighted with their 
uncertainties to obtain the sixth entry in Table [5] ("Combined F'fl 



4 The uncertainty quoted there is the dispersion of the five values, not 
the uncertainty of the mean. Even though the five values appear to be 
compatible among them (their normalized dispersion with respect to the 
weighted mean is 1.6), they are not fully independent measurements of 
the period, since the EWs and line ratios were not obtained at completely 
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Figure 6. Stokes / and V profiles in the co-added ESPaDOnS spectrum of NGC 1624-2 (phase 0.96). Left: C iv ,1,15801, 5812 lines, along with the Diffuse 
Interstellar Band (DIB) feature at 5797 A. Middle: He i ,15876. Right: He i /17066. Note the clear Stokes V signatures present in all of the stellar lines, and the 
lack of any signature in the DIB feature. Note V and N have been scaled (in the left and middle panels) and shifted for display purposes. 



Table 6. Longitudinal field measurements obtained from individual lines in the ESPaDOnS and Narval spectra of NGC 1624-2, evaluated using Eq. (2). 
Uncertainties are formal errors computed by propagation of the spectral error bars. Listed are the heliocentric Julian date of the midpoint of the observation, 
the instrument used, the peak signal-to-noise ratio per 2.6 kms~' velocity bin, and the evaluation of the longitudinal field from Stokes V. In no case is any 
detection obtained in the N profiles. The first 5 rows provide the results for the individual ESPaDOnS spectra. The 6th row provides the statistical mean and 
standard deviation for the measurements of each spectral line. The 7th row provides the longitudinal field measured from the mean spectrum obtained from 
the co-addition of the 5 individual ESPaDOnS spectra. The final row provides the results for the Narval spectrum. 



HJD 


Inst 


SNR 








Be ± 


o-s (kG) 








-2455000 




(bin- 1 ) 


C iv -15801 


Civ ,15811 


He i ,14712 


He i ,14921 


He i ,15015 


He ii ,15411 


He 1,17281 


O m ,15591 


958.715 


E 


140 


6.62 ± 1.4 


8.27 ± 1.1 


5.39 ± 1.0 


7.46 ± 1.1 


7.10± 1.3 


3.50 ± 1.1 


4.71 ±0.7 


8.91 ± 1.3 


959.716 


E 


135 


8.94 ± 1.4 


2.47 ± 1.0 


3.21 ± 1.0 


5.49 ± 1.2 


5.19 ± 1.2 


1.89 ± 1.3 


4.34 ± 0.8 


8.51 ± 1.5 


960.713 


E 


86 


4.29 ± 1.7 


5.58 ± 1.3 


3.80 ± 1.6 


9.70 ± 1.5 


9.06 ± 1.5 


1.68 ± 1.8 


6.56 ± 0.9 


7.28 ± 1.6 


961.713 


E 


137 


5.65 ± 1.4 


3.52 ± 1.0 


2.95 ± 1.0 


3.15 ± 1.1 


8.08 ± 1.1 


5.34 ± 1.2 


2.76 ± 0.8 


12.56 ± 1.4 


966.720 


E 


116 


8.47 ± 1.6 


4.06 ± 1.1 


8.92 ± 1.3 


7.84 ± 1.4 


3.96 ± 1.3 


5.85 ± 1.4 


6.67 ± 0.9 


1.63 ± 1.5 


W)±cr 


E 




6.74 ± 1.9 


4.77 ± 2.3 


4.85 ±2.5 


6.73 ± 2.5 


6.68 ±2.1 


3.65 ± 1.9 


5.00 ± 1.6 


7.78 ±4.0 


Co-added 


E 


300 


6.59 ±0.8 


4.75 ± 0.6 


4.66 ±0.7 


6.30 ± 0.6 


6.89 ± 0.7 


3.73 ±0.6 


4.73 ± 0.4 


8.98 ± 1.0 


1011.332 


N 


92 


6.50 ± 1.6 


7.66 ± 1.4 


2.67 ± 1.4 


5.78 ± 1.0 


4.80 ± 1.7 


5.51 ± 1.6 


3.67 ± 1.0 


0.80 ± 0.6 



Table 5. Periods derived from analysis of spectra line EW variations. The 
procedures used are described in the text. The adopted period is that corre- 
sponding to "Combined 1". The uncertainty quoted there is the dispersion 
of the five values, not the uncertainty of the mean. Even though the five 
values appear to be compatible amongst themselves (their normalized dis- 
persion with respect to the weighted mean is 1.6), they are not fully inde- 
pendent measurements of the period, since the EWs and line ratios were not 
obtained at completely different times for the different quantities (e.g. H/3, 
He ii ,14686, and He n ,14542 / He i ,1447 1 were obtained at exactly the same 
times). 



Quantity 
measured 



Period 
(d) 



Ha 

He i ,15876 
Hy3 

He ii ,14686 

He ii ,14542 /He U4471 
Combined 1 
Combined 2 



156.57 ±0.53 
158.38 ± 1.70 
158.97 ± 0.66 
158.52 ±0.52 
158.31 ±0.65 
157.99 ± 0.94 
156.87 ±2.95 



different times for the different quantities (e.g. the H/3, Hen ,14686, and 
He ii ,14542 / He i ,1447 1 were obtained at exactly the same times). 



Secondly, before performing the period analysis we combined the 
various datasets by subtracting the mean value from each of the five 
quantities and dividing them by their respective dispersions. Then 
we calculated the periodogram from the combined dataset, yielding 
the last entry in Table [5] ("Combined 2"). The uncertainty there 
is significantly worse than for the other five, which we attribute 
to the fact that the phased profile variations for the five measured 
quantities are not exactly the same. 

We also checked that longer periods (e.g. twice the inferred 
period, i.e. 158 x 2 = 316 d, corresponding to a double-wave vari- 
ation) were incompatible with the data (in that they generated no 
coherent variation). No other acceptable phasing is found at this 
period or any other period. 

Therefore, our preferred period is the sixth entry in Table |5] 
157.99+0.94 days, and we adopt the ephemeris of the spectroscopic 
variations as follows: 



JD = (2455967.0 ± 10) + (157.99 ± 0.94) • E. 



(1) 



According to this ephemeris, the ESPaDOnS spectropolari- 
metric observations were acquired at phases 0.95-1.0, while the 
Narval data were obtained at phase 0.28. 

If we assume that the adopted variability period corresponds 
to the period of rotation of the star (as has been confirmed through 
the o blique rotator model for other Of?p stars: IWade et aLll20 1 lL 
120121) . we tentatively conclude that NGC 1624-2 is a very slowly 
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rotating star, with a rotational period nearly 0.5 years in duration, 
and a single-wave variation of the EWs of its emission and absorp- 
tion lines. Such slow rotation would be consistent with the negli- 
gible rotational broadening of spectral lines reported in Sect.[4]and 
Sect. [6] 



6 MAGNETIC FIELD 

The individual reduced ESPaDOnS and Narval Stokes V spectra 
exhibit clear Zeeman signatures in profiles of various lines. Given 
the relatively low SNR (~ 100-150) of the individual observations, 
this implies the presence of a very strong magnetic field. 

Examination of the ESPaDOnS Stokes / spectra shows no sig- 
nificant variability of absorption or emission lines during the period 
of observation with that instrument (Feb 1-9 2012, phases 0.95-1.0 
according to Eq. (1)). This is consistent with the long (158 d) pe- 
riod of spectral variability inferred for NGC 1624-2 (Sect. |5J. Tak- 
ing advantage of the stability of the spectrum, we co-added the 5 
ESPaDOnS observations. Examples of Zeeman signatures present 
in the C iv -U5801, 5812, He i ,15876 and He i ,17065 lines of 
this mean spectrum are illustrated in Fig. [6] According to Eq. (1), 
the flux-weighted phase corresponding to the co-added ESPaDOnS 
spectrum is 0.96. 

6.1 Longitudinal magnetic field from individual lines 

We began by measuring the longitudinal magnetic field (B s ) 
of individual spectral lines in each of the polarized spec- 
tra. We selected lines in absorption that appeared to be the 
least contaminated by emission and which showed clear Stokes 
V signatures in the co-added spectrum: C iv /1/15801, 581 1, 
He i X14712, 4921, 5015, 7281, He n ,15411, and O m ,15591. We 
used the first-moment method of iRees & Semell jl979h . All lines 
in the ESPaDOnS and Narval spectra are observed to have signif- 
icant shifts from their rest wavelengths, corresponding to a radial 
velocity of about -32 kms 1 . We integrated the /// c and V/I c pro- 
files about their centres-of-gr a vity v in velocity v, in the manner 
implem ented by iDonati et al.l d 19971) and corrected by IWade et all 
fcOOOl) : 

f (v-v )V(v)dv 

<B ; > = -2.14 x 10" — . (2) 

Azc I [1 -I(v)]dv 

In Eq. (2) V(y) and I(v) are the V/I e and /// c profiles, respec- 
tively. The wavelength A is expressed in nm and the longitudinal 
field (6 ; ) is in gauss. The wavelength and Lande factor z corre- 
spond to those of each individual spectral line. Atomic data were 
obtained from the Vienna Atomic Line Database (VALD) where 
available. When experimental Lande factors were unavailable, they 
were calculated assuming L-S coupling. We integrated approxi- 
mately symmetrically about vo, using the observed span of the / 
and V profiles to establish the integration bounds. The (B-) mea- 
surements obtained from each of the lines in the individual spectra, 
and the co-added ESPaDOnS spectrum, are summarized in Table[6] 
The relevant atomic data are contained in Table [7] 

The grand mean obtained by averaging over all lines in Table|6] 
in the co-added ESPaDOnS spectrum is 5.83 kG, with a standard 
deviation of 1.7 kG. This is consistent with the longitudinal field 
obtained from averaging the measurements from all lines in the 
individual ESPaDOnS spectra (5.78 kG with a standard deviation 



of 1.4 kG). The Narval spectrum yields a similar mean value of 
(S ; >: 4.67 kG with a standard deviation of 2.2 kGThe longitudinal 
magnetic field of NGC 1624-2 is clearly remarkably strong. 

Table [6] shows that most lines exhibit somewhat more scatter 
in their measurements than is predicted by the formal errors (i.e. the 
standard deviation in row 6 is frequently ~ 2 times larger than the 
typical formal error of individual measurements). Examination of 
the measurements in greater detail revealed that the scatter results 
primarily from measurement of the first moment of V/I c , i.e. the 
numerator of Eq. (2). To explore the origin of the scatter, we first 
inspected the Stokes V profiles of a given line, which we found to 
be identical within the noise. We then investigated the sensitivity 
of the measurements to adjustments of the integration range. We fi- 
nally explored different methods of establishing vo (e.g. using only 
the Stokes V profile, using Stokes V weighted by the Stokes / line 
depth, fixing vo at a constant value). Ultimately, we concluded that 
the extra fluctuations are likely due to the difficulty of establish- 
ing the position of the centre-of-gravity v in the relatively noisy 
profiles of the individual spectra. This is supported by the mea- 
surements from the higher SNR co-added spectra (i.e. row 7 of Ta- 
ble^, which display no relatively extreme values, i.e. these values 
are in quite good agreement with the mean of the measurements 
for each line (row 6). We therefore conclude that measurements of 
(B-) from individual lines in the individual spectra are susceptible 
to relatively large fluctuations, while the averaged measurements 
(or measurements obtained from higher SNR co-added spectra) are 
more robust. 

We also observe that some spectral lines yield significantly 
lower/higher values of the longitudinal field than the mean (e.g. 
He n ,15411 (3.73 ± 0.6 kG in the co-added ESPaDOnS spectrum), 
vs. O m ,15591 (8.98 ± 1.0 kG)). This may result from several ef- 
fects. For example, most lines in the spectrum of NGC 1624-2 are 
variable. The origin of this variability is not understood in detail; it 
may be a result of changes in the formation environment of a par- 
ticular spectral line, or varying contributions of the photosphere vs. 
the magnetosphere. If lines are affected differently by these (poten- 
tially various, potentially complex) contributions, this may modify 
the numerator of Eq. (2) (e.g. if the line forms in a region where 
the magnetic field may be intrinsically weaker, e.g. the circum- 
stellar environment) or the denominator (e.g. if an absorption line 
is infilled by unpolarised emission). A thorough investigation of 
these effects will require good-quality observations spanning the 
rotational cycle of the star, coupled with a more sophisticated un- 
derstanding of the line formation processes and the origin of the 
variability. This star is uniquely well-suited to such an investiga- 
tion. 

6.2 Longitudinal magnetic field from LSD profiles 

We also analyzed the mean Stokes / and V p rofiles computed b y 
means of Least-Squares Deconvolution (LSD. lDonati et af]| 19971) . 
Donati's implementation of LSD was applied to the individual and 
co-added Stokes spectra. We took advantage of the unique charac- 
teristics of NGC 1624-2 - the detectability of Zeeman signatures 
in many individual lines - to develop a custom line mask for this 
star. The mask includes 8 spectral lines (not surprisingly, the same 
8 spectral lines analyzed above), and is summarized in Table|7] 

Using t he y 2 signal detection criteria described by 
IDonati et al.l dl997l) , we evaluated the significance of the sig- 
nal in both LSD Stokes V and in N. In no case is any signal 
detected in N, while signal in V is detected definitely (false alarm 
probability fap < 10~ 5 ) in all of our 5 ESPaDOnS observations, as 
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Figure 7. LSD profiles of NGC 1624-2. Left: LSD profiles from the first mask corresponding to absorption lines. The weak absorption feature at 
~ -200 kms~' is due to the DIB feature at 5797 A. Right: LSD profiles of weak O m emission lines. Note the inverted Stokes V profile relative to the 
left panel. This phenomenon is discussed in Sect. 6.4. 
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Figure 8. Left - Splitting of the Stokes / profiles of the C iv /U5801, 5812 lines in the Narval spectrum. Overplotted are model profiles for a ~ 20 kG dipole 
viewed at the positive magnetic pole. A solar abundance of C is assumed, and no additional broadening has been added. Right - The same theoretical profiles 
are compared to the co-added ESPaDOnS spectrum. Note the weaker lines and absence of resolved splitting in the ESPaDOnS observation. 



well as in the Narval observation. From each set of LSD profiles 
we measured the mean longitudinal magnetic field in both V and N 
using Eq. (2), integrating from -125 to +90 kms~' as determined 
from the radial velocity and extent of the Stokes V and / LSD 
profiles. We used the mean SNR-weighted Lande factor and 
wavelength, averaged over all lines in the mask. The longitudi- 
nal field measured from Stokes V is detected significantly (i.e. 
|z| = \B[\/cr > 3) in all but one of our observations, and in the mean 
ESPaDOnS spectrum at 10.7cr confidence «B ; > = 5.35 ± 0.5 kG). 
The longitudinal field in the Narval spectrum, 5.07 + 1.2 kG, 
is in formal agreement with the ESPaDOnS result obtained 
approximately 0.3 cycles earlier. In no case is the longitudinal 
field significantly detected in N. The results of the longitudinal 
field analysis are summarised in Table[8] The LSD Stokes /, V and 
N profiles extracted from the co-added ESPaDOnS spectrum are 
illustrated in Fig.|7](left frame). 

The longitudinal field inferred using LSD is formally consis- 
tent with that obtained from the individual line measurements from 
the co-added spectrum, using the standard deviation as the uncer- 
tainty. We note that the formal uncertainty associated with the LSD 
(B ; ) measurement is substantially larger than the formal error we 



would compute from combining the co-added spectrum measure- 
ments (equal to 0.2 kG). This is a consequence of the inherent scal- 
ing applied to the LSD error bars in order to obtain a satisfactory 
agreement betw een the LSD polar i zed spectrum and the real obser- 
vations (see e.g. IWade et ai]|2000l : [silvester et alj|2009r) . The LSD 
error bars therefore take into account line-to-line systematic errors 
such as those described in Sect. 6.1, and therefore probably provide 
a more realistic estimation of the uncertainties. As a consequence, 
we will adopt the LSD (B z ) measurements in the discussions that 
follow. 



The strong longitudinal field inferred from the individual lines 
and LSD profiles fully supports our initial impression that the mag- 
netic field of NGC 1624-2 is very strong. Assuming a limb dark- 
ening coefficient u = 0.3 and a centred dipolar magnetic field, 
such a strong longitudinal field implies a (minimum) polar surface 
field strength o f the dipole of nearly 20 kG according to Eq. (1) of 
IPrestonl dl967l) . Such a field is nearly 8 times stronger than that of 
any other known magnetic O-type star. 
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6.3 Zeeman splitting and magnetic field modulus 

As discussed in Sect. [4] metallic spectral lines of NGC 1624-2 are 
very sharp. Given the strong magnetic field inferred above, a signif- 
icant component of this broadening may result from the magnetic 
field. In fact, while examination of the Stokes / spectral lines in 
the ESPaDOnS spectra reveals no clear evidence of magnetic split- 
ting, the C iv /1/15801, 5812 lines in Narval spectra appear to show 
splitting in the line core (Fig. [8] left panel). 

To characterize the magnetic properties of NGC 1624-2 in 
more detail, to infer the surface field modulus (i.e. the disc- 
integrated modulus of the magnetic field, to which the Stokes / 
Zeeman splitting is most straightforwardly related) and test the re- 
ality of the large longitudinal field discussed above, we attempted 
to reproduce the Stokes / and V profiles of C iv /M5801, 581 1 using 
spectrum synthesis, taking into acco unt the influence of the mag- 
netic field. We used the Zeeman code jLandstreedf988l;IWade etall 
1200 lb . which computes profiles of spectral lines in LTE, taking 
into account Zeeman splitting and solving the equations of ra- 
diative transfer in all four Stokes parameters. We assumed a so- 
lar abundance ATLAS9 model atmosphere with T e g = 35 kK and 
log g = 4.0. The assumption of LTE is probably rather poor in the 
case of such a hot star as NGC 1624-2, and likely questionable 
for inferring detailed abundances. Nevertheless, we expect that the 
splitting properties of spectral lines are not likely to be influenced 
in any significant way. 

To begin we adopted a dipole magnetic field oriented such that 
our line-of-sight is aligned with the positive magnetic pole. Ad- 
justing the loc al polar field strength to that predicted by Eq. (1) of 
I Preston! dl967b . the Zeeman code reports a longitudinal magnetic 
field of 5.2 kG, in agreement with the measured field. The reported 
surface field modulus is 15 kG. As shown in Fig. [8] (left panel), in 
the absence of any other line broadening, this model does a rea- 
sonable job of reproducing the observed splitting and line width in 
the Narval spectrum. (We find in fact that while the best-fit field 
modulus to the /15801 line is 15 kG, for the /15812 line it is closer 
to 13 kG). This implies that the bulk of the rotational and turbu- 
lent broadening assumed in Sect.|4]is likely magnetic. The Stokes 
V profiles predicted by the model are in agreement with the ob- 
served circular polarization, although the Narval observations are 
sufficiently noisy that this agreement lacks any real meaning. From 
this comparison we conclude that the Stokes / splitting observed 
in the Narval spectrum is consistent with the surface field expected 
from a nearly 20 kG magnetic dipole, and that in the presence of 
such a field other contributions to the broadening of the C iv lines 
(e.g. rotation, turbulence) are negligible (at least at the phase of the 
Narval observation). 

Performing the same comparison with the co-added ES- 
PaDOnS spectrum (Fig. [8] right panel) we see that the C iv lines 
are weaker than in the Narval spectrum, and that no splitting is ob- 
vious. We find that this can be achieved in the models by reducing 
the field modulus and adding a small additional broadening (rota- 
tional, turbulent or magnetic) to reproduce the wings of the line 
profile. This corresponds to a total equivalent broadening of about 
30 kms~' . We find that the largest field moduli corresponding to 
dipoles that we are able to accommodate in this manner are about 
12 kG for ,15812, and 14 kG for /15801. Thus it appears that the field 
modulus inferred from the ESPaDOnS spectrum (at phase 0.96) is 
at least slightly weaker (upper limit of 13 ± 1 kG) than that inferred 
from the Narval spectrum (14 ± 1 kG, at phase 0.28). 



Table 7. Atomic data of lines used for magnetic analysis, including for con- 
struction of line masks used for Least-Squares Deconvolution. 



Wavelength 


Species 


Xlow 


Lande 


(A) 




(eV) 


factor 




Absorption 


lines 




4713.139 


He i 


20.964 


1.250 


4713.156 


He i 


20.964 


1.750 


4921.931 


He i 


21.218 


1.000 


5015.678 


He i 


20.616 


1.200 


5592.252 


On 


33.858 


1.000 


5801.312 


Civ 


37.549 


1.167 


5811.968 


Civ 


37.549 


1.333 


7281.349 


He i 


21.218 


1.200 




Emission lines 




7306.847 


Ora 


44.230 


1.200 


7307.117 


Om 


44.243 


1.750 


7455.356 


Om 


44.277 


1.200 


7515.987 


Om 


44.277 


2.250 


8172.149 


Om 


44.470 


1.000 



Table 8. Longitudinal field measurements of NGC 1624-2 obtained from 
LSD absorption line profiles. The final row provides the results for the 
mean spectrum obtained from coaddition of the 5 independent observa- 
tions. In the "Inst" column, E=ESPaDOnS, N=Narval. In the "Det?" col- 
umn, DD=Definite Detection. 



HJD 


Inst 


SNR 


Det? 


LSD B( 


LSD N( 


-2455000 




pix~' 




(kG) 


(kG) 


958.715 


E 


140 


DD 


5.95 ± 1.0 


-1.61 ±0.9 


959.716 


E 


135 


DD 


6.33 ± 1.0 


-0.48 ± 1.0 


960.713 


E 


86 


DD 


6.01 ± 1.8 


-1.30 ± 1.8 


961.713 


E 


137 


DD 


4.76 ±0.9 


-0.03 ± 1.0 


966.720 


E 


1 16 


DD 


3.34 ± 1.3 


+2.25 ± 1.3 


Co-added 


E 


300 


DD 


5.35 ±0.5 


-0.33 ± 0.5 


ESPaDOnS 












1011.332 


N 


92 


DD 


5.07 ± 1.2 


+ 1.37 ± 1.2 



6.4 Inverted Stokes V profiles of weak emission lines 

The Stokes / and V spectra of NGC 1624-2 are quite complex. For 
example, Fig.[6]shows 5 absorption/emission features. Of the three 
absorption features, one is unpolarised (the (interstellar) DIB), and 
two display Stokes V profiles consistent with a positive longitu- 
dinal field (i.e. with positive circularly-polarised flux in the blue 
wing). The two emission lines (He i /15876 and He i /17065) show 
Stokes V profiles similar to those of the absorption lines. However, 
because these lines are in emission (i.e. they have negative equiva- 
lent widths), the longitudinal field that would be inferred from these 
lines would in fact be «egaf/v<jf]- 

During our exploration of the Stokes V spectrum, we noticed 
that a number of weak emission lines of multiply-ionised light ele- 
ments, primarily O m, also exhibited Zeeman signatures. However, 
these signatures appeared to be inverted relative to those observed 

5 This complex behaviour underscored the value of high spectral resolving 
power for investigating magnetic fields of objects like NGC 1 624-2. 
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in the absorption lines. We constructed a second LSD mask consist- 
ing only of those lines (see Table |7), and extracted the mean LSD 
profile. As shown in Fig.|7](right panel), our initial suspicions were 
verified. These weak emission lines do exhibit inverted Stokes V 
signatures relative to the absorption lines and other strong emission 
lines such as He i A5876. 

According to our CMFGEN models, these lines should form 
very close to the photosphere in normal O-type stars. In the spec- 
trum of NGC 1624-2, the width and radial velocity of these lines 
are consistent with those derived from the C iv absorption lines, 
suggesting that this is probably true is the present case. 

The reversed signatures in these lines can be naturally ex- 
plained if the V profiles are formed in emission. This is in contrast 
to the other emission lines - "wind lines" - which exhibit V sig- 
natures with signs that are consistent with those of the absorption 
lines. We interpret this to indicate that these latter signatures do not 
form in emission (i.e. in the emitting region). Rather, they would 
form primarily in the photosphere in absorption lines (where the 
source function is decreasing with height). The emission contribu- 
tion of the magnetosphere and wind - dominating the observed lines 
- would presumably form at relatively large distances from the star, 
where the dipole field has decreased significantly in strength, and 
would therefore contribute only weakly to the observed V profiles. 

Measuring the longitudinal field from the O ni LSD profiles 
yields (B-) = 2.58 + 0.7 kG. The longitudinal field from these emis- 
sion lines is therefore consistent in sign with that measured from the 
absorption lines, but more than a factor of 2 smaller. However, as 
discussed above the emission lines are not predicted to probe sub- 
stantially different spatial regions from the absorption lines, at least 
as estimated from a spherical wind model. An intriguing possibil- 
ity that potentially explains the weaker longitudinal field strength 
has these lines forming primarily in the low-velocity plasma con- 
fined in closed magnetic loops, around (5.35/2.58) 1 ' 3 - 1 - 0.3 R, 
above the stellar surface. Such a scenario is qualitatively consistent 
with the expected region of formation, the measured characteristics 
of the Stokes / profiles (width and radial velocity), and the weaker 
measured longitudinal field of consistent sign. On the other hand, 
given the uncertainties related to the detailed formation of lines 
in the photosphere/magnetosphere/wind discussed in Sect. 6.1, we 
reserve further speculation until the line formation is better under- 
stood. 



7 X-RAYS 

NGC 1624-2 was observed for 10 ks with the ACIS-I CCD array on 
the Chandra X-ray Observatory in 2006 (ObsID 7473, PI Garmire, 
phase 0.12 according to Eq. 1), in which NGC 1624-2 was detected 
with 38 counts. The source is located on-axis, but very close to a 
gap between the 12 and 13 CCDs. Accordin g to the exposure map 
provided by the Chandra Source Catalogue I Evans et al.l2010l) , the 
unfortunate positioning of the source on the detector resulted in a 
loss of 35-40% in effective area. 

We extracted the source spectrum with specextract in 
CIAO 4.4, using a circular source region with a radius of 5 pixels 
(~ 2.5 arcsec), which encircles more than 95% of the energy at an 
off-axis angle of 0.09 arcmin. The background was extracted from 
a larger, source-free region nearby, straddling the same fraction of 
the two chips. Despite the poor data quality, it is nonetheless clear 
that the spectrum is very hard, with only one count below 1 keV. 
This suggests either very hot plasma, or a highly absorbed, some- 
what cooler plasma. 



Table 9. Single-temperature APEC plasma models, absorbed by a cool ISM 
model. 



N ISM 
r>22 „ m -2 



yytocal 
n22 „ m -2 



kT 



(lO^cnr 2 ) (10 2/ cm- 2 ) (keV) (1(T 14 erg s" 1 cnT 2 ) (lO^ergs" 1 ) 



Abs. flux' 

-14, 



Lx" 



0.48 
0.48 



9 ^+0.9 
z.j_ 03 

i n+u.7 n 7 +0.4 





3+0.7 



5 +2 



9 4+2.0 
1 s+1-6 



a In the 0.5- 7.0 keV band. 

* The absorbed flux was corrected for the iV"J? M only. 



We modelled the ACIS-I CCD spectrum in xspe c 12.7 iArnaudl 
1996) with a single-temperature, solar ab undance ( Asplund et al .1 



20091) APEC plasma (with atomd b 2.0.1. ISmith et alJ|200ll) . and 

t alj J2000h 



Given the 



the cold ISM absorption model of IWilmset 
low number of counts, we binned the spectrum from 0.5 to 7.0 keV 
in such a way as to obtain at least one count per bin (see Fig. [9}. 
We therefore used the C-statistic to determine the best-fit model and 
place constraints on the model parameters. We used an ISM column 
density JV{f M = 0.48xl0 22 crrT 2 , determined from E(B-V) = 0.802 
and R\j = 3.74. For such a low column density, the resulting best fit 
value for the plasma temperature is quite high (kT = 2.3^}' keV), 
although the fit is rather poor. 

We then allowed for the possibility of extra local absorption of 
the X-rays (from the wind or from a magnetosphere^), in addition 
to the ISM value of N§ M = 0.48 x 10 22 crrT 2 . The resulting fit, 
shown in Fig. [TO] is marginally better. The temperature is lower - 
kT = 0.7^2 keV - and the necessary extra column density would 
be surprisingly hi ghiV„ cal = 1.9^ x 10 22 cirr 2 compared to typical 



1_>V- JUl ^JllJlll^ljt 111 C^ll 1 1 p| X . S _(\ ^ S^. 1\J \J±±± V. V I I I pill L\_* LJjJJ.^tll. 

O-type stars (e.g. Naze et al.l201 ll) . Fig.[TOJshows the C-statistic 1, 
2 and 3 cr contours for the plasma temperature and local absorbing 
column, to illustrate the possible range in parameters. Higher SNR 
data are therefore needed to obtain robust constraints on the plasma 
temperature and absorbing column. 

The absorbed fluxes for the two models described above range 
from 3-6xl0~ 14 erg s~' cirr 2 . Correcting for the ISM column den- 
sity only and using a distance of 5.2 ± 1.1 kpc, the X-ray luminosity 
is l-4xl0 32 ergs" 1 . With a luminosity of log(L/L ) = 5.1, the X- 
ray efficiency is therefore log(Lx/i*oi) ~ -6.4, a fact or of about 
4 hig her than the canonical value for O stars of ~ -7 jNaze et all 
1201 ll) . 

High X-ray luminosity is a common feature among mag- 
netic O-type stars with confined winds. The Of?p stars HD 108, 
HD 191612 and HD 148937, as w ell as Q7 V star g' OriC, have 
log(L x /L bQ i) from -6.0 to -6.2 jChlebowskil 1 19891 ; iNaze et al.1 
l2004l : iGagne et al.ll2005l : INaze et alj|2007l , 120081) . However, these 
stars show a wide range of plasma temperatures. For the magnetic 
Of?p stars the plasma is dominated by a cool component similar to 
normal O-type stars, even if some hotter plasma is present. On the 
other hand, the X-rays of the O-type star G l Ori C are dominated 
by the hard component. At first glance, NGC 1624-2 seems to have 
hotter plasma than normal O-type stars. However, a full, detailed 
modelling of NGC 1624-2 must await higher SNR X-ray observa- 
tions. 



6 Note that the local absorbing material may have different absorbi ng prop- 
erties than assumed here (see for example [Leutenegger et al.l2O10h . 
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Figure 9. ACIS-I spectra of NGC 1624-2, binned to obtain at least one count 
per bin (black dots, error bars and shaded area). The two models described 
in Tableware shown in red (solid) and blue (dotted) lines. The C-statistics 
for the fits are 13.7 and 21.4, respectively. 
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Figure 10. C-statistic 1 , 2 and 3 <t contours for the plasma temperature (kT) 
and local absorbing column (iVg ). The best fit is indicated by a cross. 



8 THE GIANT MAGNETOSPHERE OF NGC 1624-2 

Magnetospheres of O and early B stars form through the 
channelling and confinement of an outflowing wind by the 
star's magnetic field. This magnetic control breaks the sym- 
metry of radiatively-driven winds, and therefore will influence 
observable wind diagnos tics, such as Balmer line emission 
(jLandstreet & Borra 19780 and wind resonance lines in the UV 
(Sh ore & Brownll9900 . Furthermore, material forced to flow along 
the field lines will collide near the tops of closed loops, producing a 
shock-h eated volume of plasma that will eventually cool, radiating 
X-rays dBabel & Montmerle|[T997aHbT) . Given the large mass-loss 
rate expected for a 35 kK O-type star, and the observed strong mag- 
netic field, it is highly probably that such a structure exists around 

NGC 1624-2. 

As presented bv lud- Doula & Owockl J2002I) . the global com- 
petition between the magnetic field and stellar wind can be char- 



acterized by the so-called wind magnetic confinement parameter 
r t* — ^eq^*/^ v <*" which depends on the star's equatorial sur- 
face field strength (6 ec| ), stellar radius (R*), and wind momen- 
tum (Mt'oo). For a dipolar field, one can identify an Alfven radius 
Ra - r/'^R* , representing the extent of strong magnetic confine- 
ment of the wind. Above R A , the wind dominates and stretches open 
all field lines. But below R A , the wind material is trapped by closed 
field line loops, and in the absence of significant stellar rotation 
is pulled by gravity back onto the star on a dynamical (free-fall) 
time-scale. 

To estimate the wind momentum of NGC 1624-2, we deter- 
mined the theoretical wind terminal velocity given by: 



= 2.6v csc = 2.6 



2GM(\ -T e ) 



R* 



(3) 



dLamers et al.lll995l) where T e = k c L/4kGMc is the standard Ed- 
dington parameter for electron scattering opacity in a fully ion- 
ized wind K e = 0.34 cm 2 g~'. Using the adopted parameters given 
in Table [3] we o btain v„ = 2875 kmjT 1 . We used the mass- 
loss rate recipe of IVink et ail d2000l 1200 ll) and we obtain M = 
1.6 x 10~ 7 M o yr~', roughly consistent with the other Of?p stars. 
Using a dipolar field strength of 20 kcQ, this leads to a magnetic 
confinement 77* = 1.5 x 10 4 . The closed loop region of the mag- 
netosphere therefore extends up to R A]f = 11 ± 4 R* (assuming a 
conservative 40% uncertainty in fid (i.e. ±5 kG) and a factor of 
two uncertainty in wind momentum). Therefore, even taking into 
account the uncertainties, the magnetosphere of NGC 1624-2 is in- 
trinsically large, and substantially larger than that of any other mag- 
netic O-type star (for all of which R A < 4fi*). 

In the presence of significant stellar rotation, centrifugal forces 
can support any trapped material above a Kepler co-rotation radius 
fixep = (GM/u 2 ) 1 ' 3 . This requires that the magnetic confinement 
extend beyond this Kepler radius, in which case m aterial can accu- 
mulat e to form a centrifugal magnetosphere (e.g. iTownsend et al.l 
120051) . In the case of NGC 1624-2, the slow rotational period puts 
the Kepler radius much farther out (40 R*) and no long-term accu- 
mulation of wind plasma is possible. However, the transient suspen- 
sion of circumstellar material still results in a global over-density 
in the closed loops. For O-type stars with sufficient mass-loss 
rates, the resulting dynamical magnetosphere can therefore exhibit 
strong emission in B almer recombination lines dPetit et al.l |20 111 ; 
ISundqvist et al.ll2012l) . For the huge volume of the magnetosphere 
of NGC 1624-2, we would expect partic ularly strong B a lmer e mis- 
sion. In fact, as illustrated in Fig. 2 of IWalborn et al.l d2010l) and 
Fig-QTJ NGC 1624 has - by a large margin - the strongest Balmer 
line emission of all known magnetic O-type stars. 

In the prototypical magnetic O-star G' Ori C, the hard and lu- 
mino us X-rays are se e n to b e rotationally modulated dGagne et al.l 
1 19971) . iGagne et all d2005l) used 2D magnetohydrodynamical 
(MHD) simulations to show that the high temperature plasma can 
be well explained by shocks extending up to 1.5 R*, comparable to 
its inferred Alven radius of 2R*. Therefore, the trapped wind has 
only been accelerated to a fraction of its potential terminal velocity 
before reaching the top of the closed loops. 

In a larger magnetosphere such as NGC 1624-2, the wind 
speed could in principle get higher before crossing the shock front 



7 While the maximum measured field modulus is 14 ± 1 kG, to remain con- 
sistent with 77* inferred for other magnetic O stars we employ for this pur- 
pose the equatorial field of a dipole of polar strength 20 kG inferred from 
the longitudinal field. 
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Figure 11. Peak Ha emission of NGC 1624-2 (solid black line) compared 
to peak emission of HD 191612 (dashed blue line) and 6 l Ori C (dot-dashed 
red line). 

and one could expect harder X-ray emission than in other mag- 
netic O-type stars. Furthermore, if a larger fraction of the wind gets 
shock-heated, the overall emission level could be higher. However, 
the current X-ray observations do not provide sufficient constraints 
on the plasma temperature distribution, the total X-ray flux and 
the column density to be meaningfully modelled with sophisticated 
magnetosphere models. 



9 DISCUSSION AND CONCLUSION 

We have discovered an extraordinarily strong magnetic field (maxi- 
mum mean longitudinal magnetic field (B z ) = 5.35 + 0.5 kG, corre- 
sponding to a dipole of surface polar strength ~ 20 kG) in the Of ?cp 
star NGC 1624-2 that distinguishes it qualitatively from other 
known magnetic O-type stars. Of particular interest is the presence 
of clear Zeeman signatures in individual spectral lines, and the ap- 
parent detection of resolved Zeeman splitting (corresponding to a 
maximum mean magnetic field modulus of (6) = 14 ± 1 kG). The 
detectability of Zeeman splitting should in principle allow a much 
stronger constraint on the geometry and topology of the star's mag- 
netic field as compared to any other magnetic O-type star, once 
additional observations are acquired. We also suggest that NGC 
1624-2 may be a good target for future transverse Zeeman effect 
(Stokes Q and U) observations. With its sharp, magnetically-split 
lines, the QU signatures could potentially have amplitudes compa- 
rable to Stokes V. 

Using an extensive spectroscopic data set, we performed a 
first determination of the physical properties of the star. We con- 
firm that it is a main sequence object with a mass comparable to 
those of the other Of?p stars. While the models used to infer the 
stellar properties included non-LTE effects, they did not directly 
include effects of the magnetic field. Such effects may include 
desaturation of profiles of spectral lines (modifying line blanket- 
ing), and introduction of Lorentz forces, both of which may lead 
to modification of th e hydrostatic structure of the atmosphere (e.g. 
IShulvak et al . 2009.?). While the intrinsic uncertainties associated 
with the derived quantities are probably sufficiently large to domi- 
nate the observable consequences of these phenomena, considering 
the remarkable strength of the magnetic field of NGC 1624-2 their 
potential importance should be explored in more detail. 



Indeed, the spectrum of NGC 1624-2 exhibits a number of 
peculiarities that distinguish this star from other Of?p stars. First, 
we have observed unprecedented composite profiles of the C in 
/14650 complex, with narrow and broad components. We hypoth- 
esize that the broad components a re photospheric as in normal Of 
or Ofc stars dWalborn et al.ll2010h . while the narrow components 
typical of Of?p spectra are magnetospheric. Consistent with this 
interpretation, we observe that only the narrow components appear 
to vary with phase. Moreover, the narrow C m emission does not 
disappear at minimum, so unlike most other members of this class, 
NGC 1624-2 remains Of?p at both extreme phases. We have also 
found that the absorption lines of NGC 1624-2 are very narrow, and 
that the widths of C iv /M5801, 581 1 can be reproduced essentially 
by magnetic broadening. This is in stark contrast to other Of ?p stars 
for which the profiles of these lines are clearly dominated by turbu- 
lent broadening at the level of several 10s of km s _1 . The origin of 
this difference is currently unknown. 

We have also used the strong variations of various spectral 
absorption and emission lines to infer a unique and unambiguous 
spectral variation period of 157.99 ± 0.94 d. Based on the observed 
behaviour of all other known magnetic O-type stars, it is reasonable 
to assume this to be the stellar rotation period (an assumption that 
will be tested by acquisition of additional magnetic field measure- 
ments). This implies that NGC 1624-2 rotates very slowly, spin- 
ning once in approximately one-half year. Such a conclusion is in 
good agreement with the negligible v sin inferred from modelling 
the magnetically-split line profiles. It has been established that all 
known magnetic O-type stars have diverse but relative ly long pe- 
riods of rotation, from about a week (HD 14 8937: iNaze et al.l 
l2008h to perhaps more than 50 years (HD 108; INaze et alfcOOll) . 
The most common mechanism invoked to produce such slow ro- 
tation is magnetic braking, i.e. the shedding of rotational angular 
momentum via the stellar wind and enhanced lever arm provided 
by the magnetic field. If we employ the braking model described 
bv lUd-DoulaetalJ J2009I) (Eq. 25 of that paper), we can roughly 
compute the braking timescale T spin as a function of the magnetic 
field strength, stellar mass and radius, mass loss rate and termi- 
nal velocity. Using the physical paramete rs reported in Tableland 
moment of inertia coefficient k ~ 0. 1 JClaretll2004l) . we obtain 
a spin-down time for NGC 1624-2 of 0.24 Myr. As discussed in 
Sect. 1, the estimated maximum age of the cluster NGC 1624 is 
no greater than 4 Myr, i.e. no more than ~ 17T spi „. Such an age 
is, however, easily sufficient for the star to have braked from an 
initial short rotation period to its current very long period. (For 
example, if we assume the star was initially rotating at critical 
Cent = 1-1 d), the time required to slow the rotation to 158 d would 
be just ln(158/l.l) = 4.96r spm = 1.2 Myr.) On the other hand, 
if the cluster age is significantly younger than 4 Myr this would 
place constraints on the initial rotational speed of the star (requir- 
ing a slower initial rotation) or a reconsideration of the origin of the 
current slow rotation. These conclusions are subject to the assump- 
tions and limitations of the braking model (e.g. aligned magnetic 
and rotation axes computed in 2D) and the uncertainties of the in- 
put parameters (8 d , M„, R„ M, v„ and k). 

In our analysis of the optical spectrum of NGC 1624-2, we de- 
rived a upper limit on the photospheric N abundance: ([N/H]<0.3). 
We consider this upper limit to be uncertain due to our present in- 
ability to directly include the influence of the strong magnetic field 
on the line formation in our NLTE spectrum synthesis model. Ac- 
curate knowledge of the surface N abundance represents an impor- 
tant constraint on the in terior rotation profil e of a magnetic early- 
type star. As reported bv lMevnet et al.l ( l201ll) . when magnetic brak- 
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ing occurs in a massive star characterized by internal differential 
rotation, a strong and rapid mixing occurs in layers near the sur- 
face. This results in enhanced mixing, resulting in enriched surface 
abundances of nitrogen relative to similar models with no magnetic 
braking. However, when solid-body rotation is imposed in the in- 
terior, the star is slowed so rapidly that surface enrichments are 
in fact smaller than in similar models with no magnetic braking. 
Therefore magnetic braking can enhance the surface N abundance 
or, in contrast, have little effect on it, dependin g on the internal ro- 
tation profile of the star. iMartins et al. 1 d2012l) investigated the N 
surface enrichment of 6 known magnetic O-type stars, including 3 
Of?p stars. Depending on the assumed initial rotation velocity of 
the Of?p stars, they found that that they display surface nitrogen 
abundances consistent with those of non-magnetic O stars (if their 
rotation was initially rapid, v ro( ~ 300 kms~' ), or enhanced rela- 
tive no non-magnetic O stars (if their rotation was initially modest, 
v TO t ~ a few times 10 kms~' ). Assuming NGC 1624-2 has a main 
sequence age no greater than 4 Myr, comp arison of our prelimin ary 
N abundance with the figures presented bv lMartins et al.|j2012h in- 
dicates that (1) NGC 1624-2 has a N abundance slightly lower than 
non-magnetic O stars with similar positions on the HR diagram; 
(2) that the abundance is low relative to that expected for a star of 
this age and mass if it was rotating initially at high velocity; (3) but 
consistent with the expected abundance if the star was rotating ini- 
tially at lower velocity. However, as discussed above these results 
are quite tentative, and require more detailed modelling and a more 
robust determination of the age of the star before firm conclusions 
can be drawn. 

As a consequence of its intense magnetic field, NGC 1624-2 is 
expected to host a magnetospheric volume substantially larger than 
any other magnetic O-type star. The inferred magnetic wind con- 
finement parameter, 77* = 1.5 x 10 4 , is 300 times larger than that of 
the Of?p star with the next-strongest field, HD 191612. This leads 
to a predicted Alfven radius of 1 1.4 R+ (versus 2.2 R, in the case of 
HD 191612). This much larger volume of confined plasma should 
result in much stronger magne tospheric emission (e. g. according 
to the mechanisms discussed bv lSundqvist et al.l2012h . Indeed, the 
Hey emission of NGC 1624-2 is found to be substantially stronger 
than that observed in any other magnetic O-type star. For compar- 
ison, and as illustrated in Fig.QT| the ma ximum EW of Ho - in the 
spectrum of HD 191612 is a bout 4 A (eg. I Wade et al.l201 lh . while 
that of 1 Ori C is about 2 A JStahl et al.l2008l) . The peak EW of the 
Ha line of NGC 1624-2 is 26 A, 6.5 times greater than HD 191612 
and 13 times greater than that of 9 [ Ori C. Modeling of this re- 
markable Ha emission is an urgent priority, and first attempts are 
underway by Sundqvist, ud Doula et al. (priv. comm.). 

It is expected that such strong magnetic wind confinement 
in the presence of such a powerful wind should lead to intense 
X-ray emitting shocks. Analysis of archival Chandra ACIS-I X- 
ray observations indicates a hard and luminous X-ray spectrum 
(log(Lx/Lboi) 6.4), qualitatively consistent with theoretical ex- 
pectations as well as the behaviour of 9 [ Ori C. However, the cur- 
rent observations are more or less equally consistent with a rela- 
tively hard, weakly extinguished source, or a relatively soft, highly 
extinguished source. New higher-quality X-ray observations will 
be required to draw useful quantitative conclusions about the X-ray 
properties of NGC 1624-2. 

Due to the paucity of magnetic data, the magnetic topology 
and geometry of NGC 1624-2 are at present only weakly con- 
strained. Although the geometry (i.e. the inclination angle i and 
the magnetic obliquity ji) is normally inferred from the variation of 
the longitudinal magnetic field, observations of other magnetic O 



stars (e.g. 6 l Ori C, HD 191612, HD 57682) show that a clear rela- 
tionship exists between the EW variations of optical emission lines 
diagnostic of the wind (e.g. Ho-, He 1 /15876, He n /14686) and the 
longitudinal field variation. In particular, the emission extrema of 
these stars all correspond to extrema of the (sinusoidally- varying) 
longitudinal field, with the emission maximum corresponding to 
the maximum unsigned longitudinal field. If we assume that the 
magnetic topology of NGC 1624-2 is roughly dipolar, we can ex- 
pect with reasonable confidence that the ESPaDOnS observation, 
acquired at phase 0.96 near emission maximum, corresponds to 
the approximate maximum of the longitudinal field. The longi- 
tudinal field measured from the Narval observation, acquired at 
phase 0.28, is not very different from that of the ESPaDOnS mea- 
surement, suggesting that the variation of the longitudinal field is 
not very large, and moreover that it does not change sign as the 
star rotates. This implies that we view essentially only one mag- 
netic hemisphere during the rotation of the star, and consequently 
/ + y6 < 90°. This conclusion is supported by the single-wave na- 
ture of the EW variations. When both magnetic hemispheres are 
visible duri ng the stellar rotation (i.e. / + /? » 90°, as in the case of 
HD 57682; iGrunhut et alj2009l , and MNRAS, submitted), the EW 
variations of the emission lines exhibit a double-wave variation. In 
contrast, smaller values of the su m + /3 produce single-wave vari- 
ation s (e.g. 8> Ori C, HD 191612; IStahl et al.ll2008l ; lHowarth et al.l 
l2007h . 

Clearly, NGC 1624-2 holds a special place amongst the known 
magnetic O-type stars. Its extreme magnetic field and puzzling 
spectral peculiarities have much to teach us, and demand imme- 
diate attention. Plans for observational follow-up - including op- 
tical spectroscopy, spectropolarimetry, and photometry; UV spec- 
troscopy; and X-ray spectroscopy - are already underway, as are 
theoretical investigations of its Ha and X-ray properties to better 
understand the magnetospheric geometry, structure and energetics. 
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